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Abstract 



C^ , This text is a survey of recent results obtained by the author and 

cohaborators on different problems for non-self-adjoint operators. The 
topics are: Kramers-Fokker-Planck type operators, spectral asymp- 
totics in two dimensions and Weyl asymptotics for the eigenvalues of 
non-self-adjoint operators with small random perturbations. In the 
introduction we also review the notion of pseudo-spectrum and its 
f^ . relation to non-self-adjoint spectral problems. 
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1 Introduction 

1.1 Some background 

For self-adjoint and more generally normal operators on some complex Hilbert 
space Ti we have a nice theory, including the spectral theorem and a nice es- 
timate on the norm of the resolvent: 

||(z-P)-^|| < (dist(2,a(P)))~\ (r{P)= the spectrum of P. (1.1) 

This has a consequence for the corresponding evolution problem: If cr(P) C 
{z eC]^z> Ao}, then 

||e~*^|| < e^^«*, t > 0. (1.2) 

However, non-normal operators appear frequently: Scattering poles. Convection- 
diffusion problems, Kramers- Fokker- Planck equation, damped wave equa- 
tions, linearized operators in fluid dynamics. Then typically, \\{z — P)"^|| 
may be very large even when z is far from the spectrum. This implies math- 
ematical difficulties: 

- When studying the distribution of eigenvalues, 

- When studying functions of the operator, like e~*^ and its norm. 

It also implies numerical difficulties like: 

- Eigenvalue instability. 

There are (in the author's opinion) two ways out: 

- Change the Hilbert space norm to make the operators look more normal. 
(Complex scaling methods.) 

- Recognize that the region of the z-plane where ||(2; — P)~^|| is large, has 
its own interest. (Pseudospectrum.) 

The option to choose depends on the problem. 



• In some problems, like those related to scattering poles, there is no 
obvious choice of Hilbert space and we are free to make the most natural 
one. This option is particularly natural when considering a differential 
operator with analytic coefficients. 

• In other problems the canonical Hilbert space is L^ and we are at most 
allowed to change the norm into an equivalent one. Here the notion of 
pseudospectrum is likely to be important. 

Let P : H ^ H he closed, densely defined, H a complex Hilbert space 
and let p{P) = C\a{P) denote the resolvent set. The notion of pseudospec- 
trum is important in numerical analysis and we refer to L.N. Trefethen |52] . 
Trefethen-M. Embree [51] and further references given there. Thanks to 
works of E.B. Davies [5], [7], M. Zworski [57] and others it has become pop- 
ular in the non-self-adjoint spectral theory of differential operators. 

Definition 1.1 Let e > 0. The e-pseudospectrum of P is 

a,(P) := a{P) U {z E p(P); \\{z - P)-'\\ > 1/e}. 

Unlike the spectrum, the pseudospectrum will in general change when we 
change the norm on 7i. Moreover, it can be characterized as a set of spec- 
tral instability as follows from the following version of a theorem of Roch- 
Silberman 



Theorem 1.2 



a,(P)= U aiP + Q). 



QeC{H,H) 



In his survey ^52j L.N. Trefethen discusses some linearized operators from 
fluid dynamics: 

• Orr-Sommerfeld equation (Orzag, Reddy, Schmid, Hennigson). 

• Plane Poiseuille flow (L.N and A.N Trefethen, Schmid). 

• Pipe Poiseuille flow (L.N and A.N Trefethen, Reddy, Driscoll), 

and to what extent stability can be predicted from the sudy of the spec- 
trum of these non-self-adjoint operators: Eigenvalue analysis alone leads in 
some cases to the prediction of stability for Reynolds numbers R < 5772. 
Experimentally however, we have stability only for R < 1000. 

The rough explanation of this is that the e-pseudospectrum (for a suitable 
e) crosses the imaginary axis before the spectrum does, when R increases. 



Then ||e^*''^|| will grow fast for a limited time even though the growth for 
very large times is determined by the spectrum. However, since P appears 
as a linearization of a non-linear problem, that suffices to cause instability. 

In the case of differential operators the pseudospectral phenomenon is 
very general and related to classical works in PDE on local solvability and 
non-hypoellipticity. E.B. Davies [3] studied the non-self-adjoint semiclass- 
cial Schrodinger operator with a smooth (complex-valued) potential in di- 
mension 1 and showed under "generic" assumptions that one can construct 
quasimodes with the spectral parameter varying in an open complex set, con- 
taining points that are possibly very far from the spectrum (as can be verified 
in the case of the complex harmonic oscillator). M. Zworski [57] observed 
that this is essentially a rediscovery of an old result of Hormander [321 E3] , 
and was able to generalize considerably Davies' result by adapting the one of 
Hormander to the semi-classical case. With N. Dencker and M. Zworski [9] 
we also gave a direct proof and a corresponding adaptation of old results of 
Sato-Kawai-Kashiwara [44j to the analytic case: 



Theorem 1.3 (l5^, J^) Let 

dx 



P{x,hD,) = Y, ao.{x){hD,r, D, = ^ (1.3) 



|«|<-m 



have smooth coefficients in the open setVt C R". Putp{x,$,) = X]|a|<m'^a(^)^" 
Assume z = p(xo,^o) with the Poisson bracket ]^{p,p}{xo,^o) > 0. Then 3 
u = Uhe C^in), with \\u\\ = I, II (P - z)u\\ = 6{h^), when h^O. 
Analytic case: Can replace "h°°" by "e~^/'-"^". 

Here, we have used standard multi- index notation: ^" = ^"^ ■ ... ■ ^"", 
\a\ = ai + ... + an, and the norms || ■ || are the ones of L^ or i"^ if nothing 
else is indicated. 

This result was subsequently generalized by K. Pravda-Starov, [12]. No- 
tice that this implies that when the theorem applies and if the resolvent 
(P — z)~^ exists then its norm is greater than any negative power of h when 
h^O. 

Example 1.4 Let P = -h^A + V{x), p{x,^) = ^'^ + V{x), ^'^ = ^1 + ...+^l 
Then ^{p,p} = -A^ ■ '^V'{x). 

"Generically" , if 2; = p{x,^), then {p,p}{x,$,) 7^ and one can show 
quite generally that if this happens then there is also another point {y, rf) 
with p{y,ri) = z such that {p,p}{x,^) takes the opposite sign. This justifies 
the following simplified terminology in the semi-classical limit: 



The semi-classical pseudospectrum of P is the range TZ{p) of p. 

In [H] we also showed under suitable assumptions (inspired from scat- 
tering theory and from the theory of sub-elliptic operators), that there is 
no spectrum near the boundary of the semi-classical pseudospectrum and 
that we may have quite a good control of the norm of the resolvent there. 
Generalizations to the case of systems were given by Dencker [S] . 

1.2 The topics of this survey 

We will discuss three subjects involving non-self-adjoint differential and pseu- 
do differential operators, We will always wok in the semi-classical limit, which 
means that our operators are of the form P{x, hD^; h), where P is a suitable 
symbol and < /i ^ 1. It is quite clear however that some of our results will 
also apply to non-semiclassical situations in the limit of large eigenvalues. 
The subjects are: 

• The Kramers-Fokker-Planck operator, 

• Bohr-Sommerfeld rules in dimension 2, 

• Weyl asymptotics for non-self-adjoint operators with small random per- 
turbations. 

and most of the works discussed are the results of collaborations with A. Melin, 
M. Hitrik, F. Herau, C.Stolk, S. Vu Ngoc and M. Hager. 

In the first two topics we exploit the possibility of changing the Hilbert 
space norm by introducing exponential weights on phase space. 

In the case of the Kramers-Fokker-Planck operator, we make no analytic- 
ity assumptions and the phase-space weights are correspondingly quite weak. 
In this case however our operator is a differential one, so we are allowed to ap- 
ply strong exponential weights depending only on the base variables, and this 
is important when studying small exponential corrections of the eigenvalues 
via the so called tunnel effect. 

For the Bohr-Sommerfeld rules, we make analyticity assumptions that 
allow stronger phase space weights. In both cases the effect of the exponential 
weights is to make the operator under consideration more normal. 

In the third topic, we do not use any deformations of the given Hilbert 
space, but exponential weights play an important role at another level, 
namely to count zeros of holomorphic functions with exponential growth. 

The pseudospectrum will not be discussed explicitly below. In the Kramers- 
Fokker-Planck case, the problems are located near the boundary of the semi- 
classical pseudospectrum, and it turns out that we have a very nice control 



of the resolvent there. In the 2 dimensional Bohr-Sommerfeld rules, we have 
stronger exponential weights, reflecting stronger pseudospectral phenomena. 
Finally in the subject of Weyl asymptotics, we often have strong pseudospec- 
tral behaviour for the unperturbed operator. From the proofs it appears that 
the random perturbations will weaken the pseudospectral behaviour and this 
might have very interesting consequences for the associated evolution prob- 
lems. This is still very much an open problem. 

2 Kramers- Fokker-Planck type operators, spec- 
trum and return to equilibrium 

2.1 Introduction 

There has been a renewed interest in the problem of "return to equilibrium" 
for various 2nd order operators. One example is the Kramers-Fokker-Planck 
operator: 

P = yhd,- V\x) ■ hdy + ^{-hdy + y) ■ {hdy + y), (2.1) 

where x, y G R" correspond respectively to position and speed of the particles 
and h > corresponds to temperature. The constant 7 > is the friction. 
(Since we will only discuss L^ aspects we here present right away an adapted 
version of the operator, obtained after conjugation by a Maxwellian factor.) 
The associated evolution equation is: 

{hdt + P)u{t,x,y) = 0. 

Problem of return to equilibrium: Study the rate of convergence of u{t,x,y) 
to a multiple of the "ground state" uo{x, y) = e~^^ /2+v{x))/h ^j-^gj^ ^ __^ +00, 
assuming that V{x) — > -f-cxo sufficiently fast when a; ^ cxd so that mq € 
L^(R^"'). Notice here that P{uo) = and that the vector field part of P is 
h times the Hamilton field of y^/2 -|- V{x), when we identify R^^^ with the 
cotangent space of R". 

A closely related problem is to study the difference between the first 
eigenvalue (0) and the next one, /i(/i). (Since our operator is non-self-adjoint, 
this is only a very approximate formulation however.) 

Some contributions: L. Desvillettes-C Villani [10], J. P. Eckmann-M. Hairer 
, F. Herau-F. Nier [2Z!, B. Helffer-F. Nier ^2], Villani |SS|. In the work 
precise estimates on the exponential rates of return to equilibrium were 
obtained with methods close to those used in hypoellipticity studies and this 
work was our starting point. With Herau and C.Stolk ^28j we made a study in 
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the semi-classical limit and studied small eigenvalues modulo 0{h°°). More 
recently with Herau and M. Hitrik [23] we have made a precise study of 
the exponential decay of ^{h) when V has two local minima (and in that 
case /i(/i) turns out to be real). This involves tunneling, i.e. the study of 
the exponential decay of eigenfunctions. As an application we have a pre- 
cise result on the return to equilibrium [26j . This has many similarities with 
older work on the tunnel effect for Schrodinger operators in the semi-classical 
limit by B. Helffer-Sjostrand [23112^ and B. Simon ^5] but for the Kramers- 
Fokker-Planck operator the problem is richer and more difficult since P is 
neither elliptic nor self-adjoint. We have used a supersymmetry observation 
of J.M. Bismut p] and J. Tailleur-S. Tanase-Nicola-J. Kurchan [51j, allowing 
arguments similar to those for the standard Witten complex 



2.2 Statement of the main results 

Let P be given by i^^ where V G C°°(R"; R), and 

9"V(x) = C(1), |a| >2, (2.2) 

\VV{x)\ >1/C, \x\ >C, (2.3) 

l^ is a Morse function. (2.4) 

We also let P denote the graph closure of P from iS(R^") which coincides with 
the maximal extension of P in L^ (see [271 [22l [26] ) . We have 3fJP > and 
the spectrum of P is contained in the right half plane. In [28] the spectrum 
in any strip < ^z < Ch (and actually in a larger parabolic neighborhood 
of the imaginary axis, in the spirit of [27]) was determined asymptotically 
mod ((9(/i°°)). It is discrete and contained in a sector \'^z\ < C^z + 0{h°°): 

Theorem 2.1 The eigenvalues in the strip < ^z < Ch are of the form 

A,,,(/i) ~ h{fi,,k + h^'^'^^'Hk,! + h'"''-'^li,,k,2 + ..) (2.5) 

where ^j^k ore the eigenvalues of the quadratic approximation ( "non- self adjoint 
oscillator" ) 



y-d^- V"{xj)x ■ dy + -{-dy + y) ■ {dy + y), 



2 



at the points (xj, 0), where Xj are the critical points of V. 



The Hj^k are known explicitly and it follows that when Xj is not a local 
minimum, then ^Xj^k > h/C for some C > 0. When Xj is a local minimum, 
then precisely one of the Xj^k is 0{h°°) while the others have real part > h/C . 
Furthermore, when V -^ +00 as x ^ 00, then is a simple eigenvalue. In 
particular, if V has only one local minimum, then 

inf 3fJ(a(P) \ {0}) ~ /i(/ii + h^i2 + ...), /ii > 0. 

(or possibly an expansion in fractional powers) and we obtained a corre- 
sponding result for the problem of return to equilibrium. It should be added 
that when ^j^ is a simple eigenvalue of the quadratic approximation then 
Nj^k = 1 so there are no fractional powers of h in (12. 5p . 

The following is the main new result that we obtained with F. Herau and 
M. Hitrik in i25i: 



Theorem 2.2 Assume that V has precisely 3 critical points; 2 local minima, 
x±i and one "saddle point", xo of index 1. Then for C > sufficiently large 
and h sufficiently small, P has precisely 2 eigenvalues in the strip < '^z < 
h/C , namely and fi{h), where ^{h) is real and of the form 

/i(/i) = h{ai{h)e~^^'/'' + a_i(/i)e-2^-^/'^), (2.6) 

where Qj are real, 

aj{h) ~ ajfi + haj^i + ..., /i ^ 0, Qj^q > 0, 

S, = V{xo)-V{x,). 

As for the problem of return to equilibrium, we obtained the following 
result with F. Herau and M. Hitrik in 



Theorem 2.3 We make the same assumptions as in Theorem \2.2\ and let Uj 
be the spectral projection associated with the eigenvalue fij, j = 0, 1, where 
fiQ = 0, fii = fi{h). Then we have 

Uj = 0(1) ■.L'^^L\ h^ 0. (2.7) 

We have furthermore, uniformly as t > and h ^ 0, 

e-^PI^ = Ho + e-*'^^/'^Hi + 0(1)6"*/^, m L{1? , L^), (2.8) 

where C > is a constant. 

Actually, as we shall see in the outline of the proofs, these results (as well 
as (12.51) ) hold for more general classes of supersymmetric operators. 
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2.3 A partial generalization of |28] 
Consider on R" (2n is now replaced by n): 

- ^{cj{x)hdx^ + hdx^ o Cj{x)) + po{x) 
j 
= P2 + 1P1 + P0, 

where bj^k,Cj,Po are real and smooth. The associated symbols are: 
P{x,0 = P2{x,0+m{x^O+Po{x), 



P2 = J2^^'''^^^''^ P^ = J2^^^r 



Assume, 



P2 > 0, po > 0, 




d^b,,, = 0{1), \a 


>o 


d^c, = 0{1), \a 


> 1 


dy, = 0(1), \a 


>2 



Assume that 

{x; po{x) = ci{x) = .. = Cn{x) = 0} 

is finite = {xi, ...,xn} and put C = {pi, ..., p„}, pj = {xj, 0). Put 



Pix,0 = iO'W^.O+Poix), (0 = V^^+W 

1 f'^"^'^ ^ 

{p)to = 7^ P° exp{tHp^)dt, To > fixed. 

^0 J -To/2 

Here in general we let iJ^ = a^ ■ ^ — a^ ■ ^ denote the Hamilton field of the 
C^-function a = a{x,^). 

Dynamical assumptions: Near each pj we have {p)to ~ \p ~ Pjl"^ ^^id in 
any compact set disjoint from C we have {p)to ^ 1/C. (Near infinity this last 
assumption has to be modified slightly and we refer to [23] for the details.) 
The following result from j25j is very close to the main result of [2H] and 
generalizes Theorem 12. 1[ 

Theorem 2.4 Under the above assumptions, the spectrum of P is discrete 
in any band < ^z < Ch and the eigenvalues have asymptotic expansions 
as in ^2.5\) . 
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Put 

Qix, = -Pi^^ ^0 = P2ix, + Pi (a;, - Poix). 

The linearization of the Hamilton field Hg at pj (for any fixed j) has 
eigenvalues ±ak, k = l,..,n with real part ^ 0. Let A+ = A+j be the 
unstable manifold through pj for the ifg-fiow. Then A+ is Lagrangian and 
of the form ^ = (f)'^{x) near Xj (0+ = <f)+,j), where 



[X 



3/ 



0, 0'+(x,) = 0, (t>l{x,) > 0. 



The next result is from |25j : 

Theorem 2.5 Let Xj^k{h) be a simple eigenvalue as in Ii2.5\) and assume 
there is no other eigenvalue in a disc D{Xjk, h/C) for some C > 0. Then, in 
the L^ sense, the corresponding eigenfunction is of the form e~'^+^^^^'^{a{x; h) + 
0{h°°)) near Xj, where a{x; h) is smooth in x with an asymptotic expansion 
in powers of h. Away from a small neighborhood of Xj it is exponentially 
decreasing. 

The proof of the first theorem uses microlocal weak exponential estimates, 
while the one of the last theorem also uses local exponential estimates. 

2.4 Averaging and exponential weights. 

The basic idea of the proof of Theorem 12 .41 is taken from [28], but we reworked 
it in order to allow for non-hypoelliptic operators. We will introduce a weight 
on T*!^ of the form 



^e 



j J{^^)p,oexpitH,,)dt, (2.9) 



for < e ^ 1. Here J{t) is the odd function given by 

and we choose Pe{p) to be equal to p{p) when dist (p, C) < e, and fiatten out 
to ep away from a fixed neighborhood of C in such a way that Pe = 0{e). 
Then 

^pi^.= (Pe)To-Pe. (2.11) 

We let e = Ah where A ^ 1 is independent of h. Then the weight 
exp{;ip^/h) is uniformly bounded when h ^ 0. Indeed, ipt = 0{h). 

Using Fourier integral operators with complex phase, we can define a 
Hilbert space of functions that are "microlocally (9(exp(?/'e//i)) in the L^ 
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sense". The norm is uniformly equivalent to the one of L^, but the natural 
leading symbol of P, acting in the new space, becomes 

p{exp{zH^J{p)), pGT*R" (2.12) 

which by Taylor expansion has real part ~ P2{p) + Po{p) + (Pe) — Pe- 

Very roughly, the real part of the new symbol is > e away from C and 
behaves like dist(p, C)^ in a ^e-neighborhood of C. This can be used to 
show that the spectrum of P (viewed as an operator on the weighted space) 
in a band < ^z < e/C comes from an y^-neighborhood of C. In such a 
neighborhood, we can treat P as an elliptic operator and the spectrum is 
to leading order determined by the quadratic approximation of the dilated 
symbol (I2.12p . This gives Theorem 12.41 

We next turn to the proof of Theorem 12.51 and we work near a point 
Pj = {xj,C,j) G C. Recall that A+ : ^ = (p'^ix) is the unstable manifold for 
the Hg-fiow, where q{x,^) = —p{x,i^). We have q{x,(j)'_^_{x)) = 0. 

In general, if ■?/' G C°° is real, then P^ := e"^^^ o P o e~'^^^ has the symbol 

Pv(^' = P2{.x, C)-q{x, il)\x)) + i(g^(x, ip\x)) ■ i (2.13) 

• As long as q{x,ip'{x)) < 0, we have ^p^ > and we may hope to 
establish good apriori estimates for P^. 

• This is the case for ip = and for tp = 0+. Using the convexity of q{x, ■), 
we get suitable weights ip with q{x,ip'{x)) < 0, equal to 4>+{x) near Xj, 
strictly positive away from Xj and constant outside a neighborhood of 
that point. 

• It follows that the eigenfunction in Theorem l2.5l is (roughly) 0{e~'^+^^'^^'^) 
near Xj in the L^ sense. 

• On the other hand, we have quasi-modes of the form a{x; h)e~'^+^^'^^^ 
as in 



Applying the exponentially weighted estimates, indicated above, to the 
difference of the eigenfunction and the quasi-mode, we then get Theo- 
rem [2]5l 



2.5 Supersymmetry and the proof of Theorem 12.2 



We review the supersymmetry from [T], [HI], see also G. Lebeau ^36j. Let 
A{x) : T^R" — *> Tj.R" be linear, invertible and smooth in x. Then we have 
the nondegenerate bilinear form 

{u\v)a{x) = {A^A{x)u\v), u,ve a''T*R"\ 
11 



and we also write {u\v)a{x) = {u\v)a{x)- 

If u, V are smooth fc-forms with compact support, put 



{u\v)a= I {u{x)\v{x))A(x)dx. 

The formal "adjoint" Q^'* of an operator Q is then given by 

{Qu\v)a = {u\Q^^*v)a. 

Let (f) : R" — > R be a smooth Morse function with 9"0 bounded for 
|a| > 2 and with |V0| > 1/C for |a;| > C. Introduce the Witten-De Rham 
complex: 

d^ = e~h ohdoeh = "^{hd^^ + d^^cp) o dx^j 






where d denotes exterior differentiation and dx^ left exterior multiplication 
with dxj. The corresponding Laplacian is then: — A^ = d^'*d^ + d^d,'*. Its 
restriction to g-forms will be denoted by — Aj^ . Notice that: 

-A5)(e-^/^) = 0. 

Write A = B + C with 5* = B, C^ = -C . -A^ is a second order 
differential operator with scalar principal symbol in the semi-classical sense 
(fa!- ^0) of the form: 

j,k j,k 

Example. Replace n by 2n, x by {x,y), let 

1 / / 



Then 

-A J) = M0;-5x-0L-5.) 

i 

When = y^/2 + l^(x) we recover the KFP operator fl2.ll) 

The results of Subsection 12.31 apply, if we make the additional dynamical 
assumptions there; — A^ has an asymptotic eigenvalue = o{h) associated to 
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the critical point Xj precisely when the index of Xj is equal to q (as for the 
Witten complex and analogous complexes in several complex variables). In 
order to cover the cases g > we also assume that 

A = Const. (2.14) 

The Double well case. Keep the assumption (I2.14p . Assume that </> is a Morse 
function with |V0| > 1/C for |x| > C such that — A^ satisfies the extra 
dynamical conditions of Subsection 12.31 and having precisely three critical 
points, two local minima U±i and a saddle point Uq of index 1. 

Then — A)^ has precisely 2 eigenvalues: 0, ft that are o{h) while — A)^ has 
precisely one such eigenvalue: fi. (Here we use as in the study of the Witten 
complex, that d^ and d,'* intertwine our Laplacians in degeree and 1. The 
detailed justification is more complicated however.) e''^'^ is the eigenfunction 
of A^ corresponding to the eigenvalue 0. Let Sj = 4>{Uq) — (piUj), j = 
±1, and let Dj be the connected component of {x G R"; 0(x) < 4>{Uo)} 
containing Uj in its interior. 

Let E^"^^ be the corresponding spectral subspaces so that diiaE^'^^ = 2, 
dimi?*^^) = 1. Truncated versions of the function e"'^*-^^/'^ can be used as 
approximate eigenfunctions, and we can show: 

Proposition 2.6 -E''-'^^ has a basis ei,e_i, where 

ej = xj(x)e-^('^(^)-'^(^^» + 0{e~r^(S'-'^), m the L^-sense. 
Here, we let Xj G C^{Dj) be equal to 1 on {x E Dj] (f){x) < (piUo) — e}. 
The theorems 12.41 12.51 can be adapted to — AJ4 and lead to: 
Proposition 2.7 E^^"* = Ceo, where 

eo{x) = Xoix)aoix; /i)e-^'^+(^) + 0(e-^«/'^), 

0+(x) ~ (x — f/o)^, Co > is small enough, oq is an elliptic symbol, xo £ 
Co-(R"), Xo = l nearUo. 

Let the matrices of (i^ : E^^'^ -^ E^^^ and d^'* : E^^^ — > E^^'^ with respect 
to the bases {e_i,ei} and {cq} be 

( A_i Ai ) and I ,^^ J respectively. 

Using the preceding two results in the spirit of tunneling estimates and com- 
putations of Helffer-Sjostrand ([231 [21]) we can show: 
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Proposition 2.8 Put Sj = (p{Uo) — (p{Uj), j = ±1. Then we have 

^-l\ z^i/r ^^/ -^^^ f i-lih)e-^-^/^ 

X,')=h'i' + ^i^ ^^))( il\h)e~^./^ 

^*l\ z^i/r ^^/ -^^^ f tAh)e-^-''^ 

where i±i, i^i are real elliptic symbols of order such that iji* > 0, j = ±1. 



From this we get Theorem I2.3[ since jj. = \*_^\_i + A^Ai. D 

Thanks to the fact that we have only two local minima, certain simpli- 
fications were possible in the proof. In particular it was sufficent to control 
the exponential decay of general eigenfunctions in some small neighborhood 
of the critical points. For more general configurations, it might be necessary 
to get such a control also further away and this seems to lead to interesting 
questions, involving degenerate and non-symmetric Finsler distances. 

2.6 Return to equilibrium, ideas of the proof of Theo- 
rem 12.3 



Keeping the same assumptions, let Ho, Hi be the rank 1 spectral projections 
corresponding to the eigenvalues /io := 0, /ii := /i of — A^ and put 11 = 
Ho + Hi. Then e_i,,ei is a basis for 7?.(n) and the restriction of P to this 
range, has the matrix 

'^; ) ( ^-. ^. ) - ( "^t: %': ) (2-1=) 

with the eigenvalues and /i = Al^A-i + A^Ai. A corresponding basis of 
eigenvectors is given by 

vo = ^(Aie_i-A_ie_i) (2.16) 



vi = ^=(Alie_i + A*e_i). 



The corresponding dual basis of eigenfunctions of P* is given by 

V* = l-{XleU - y.eU) (2.17) 

vl = ^(A_ie*_i + Aie*), 
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where e*_i,e\ G 7^(11*) is the basis that is dual to e_i,ei. It follows that 
Vj, V* = C(l) in L^, when h ^ 0. 

From this discussion we conclude that IIj = {■\v*)vj, are uniformly bounded 
when h ^ 0. A non-trivial fact, based on the analysis described in Subsec- 
tions 12.31 12.41 is that after replacing the standard norm and scalar product 
on L^ by certain uniformly equivalent ones, we have 

^{Pu\u)>-\\uf, VMG7^(l-^), (2.18) 

where 11 is the spectral projection corresponding to the spectrum of P in 
D{0, Bh) for some S > 1. 

This can be applied to the study of u{t) := e^*^/^-u(0), where the initial 
state u(0) G L^ is arbitrary: Write 

m(o) = noM(o) + niM(o) + (1 - n)M(o) =: m° + m^ + u^. (2.19) 

Then 

h°IUki,||«^|| < O{l)\\u{0)\\ (2.20) 

||e~*^/\^|| < Ce-*/^||M(0)|| (2.21) 

e-*^/%- = e-*^^%-, J = 0, 1. (2.22) 

Here (E^Il) follows if we write m^ = (1 - n)^ + (H - Il)u, apply fl2A8D to 
the evolution of the first term, and use that the last term is the (bounded) 
spectral projection of it to a finite dimensional spectral subspace of P, for 
which the corresponding eigenvalues all have real part > h/C. □ 

3 Spectral asymptotics in 2 dimensions 

3.1 Introduction 

This section is mainly based on recent joint works with S. Via Ngoc and 
M. Hitrik [22] [30] , but we shall start by recalling some earlier results that we 
obtained with A. Melin [10] where we discovered that in the two dimensional 
case one often can have Bohr-Sommerfeld conditions to determine all the 
individual eigenvalues in some region of the spectral plane, provided that 
we have analyticity. This was first a surprise for us since in the self-adjoint 
case such results are known only in 1 dimension and in very special cases for 
higher dimensions. 

Subsequently, with M. Hitrik we have studied small perturbations of self- 
adjoint operators. First we studied the case when the classical flow of the 
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unperturbed operator is periodic, then also with S. Vii Ngoc we looked at the 
more general case when it is completely integrable, or just when the energy 
surface contains some invariant diophantine Lagrangian tori. 

3.2 Bohr-Sommerfeld rules in two dimensions 

For (pseudo-) differential operators in dimension 1, we often have a Bohr- 
Sommerfeld rule to determine the asymptotic behaviour of the eigenvalues. 
Consider for instance the semi-classical Schrodinger operator 

P = -h^-r^ + Vix), with symbol p(x, ^) = ^'^ + Vix), 

where we assume that V G C°°(R; R) and V{x) -^ +cxd, |x| -^ cxd. Let Eq G 
R be a non-critical value of V such that (for simplicity) {x G R; V{x) < Eq} 
is an interval. Then in some small fixed neighborhood of Eq and for h > 
small enough, the eigenvalues of P are of the form E = Ek, k & Z, where 



k-9{E;h), I{E) = f ^ ■ dx, 9{E;h) r^ 9o{E) + 9i{E)h 

Jp-HE) 



2Txh Jp-^{E) 

In the non-self-adjoint case we get the same results, provided that ^V is 
small and V is analytic. The eigenvalues will then be on a curve close to the 
real axis. 

For self-adjoint operators in dimension > 2 it is generally admitted that 
Bohr-Sommerfeld rules do not give all eigenvalues in any fixed domain except 
in certain (completely integrable) cases. Using the KAM theorem one can 
sometimes describe some fraction of the eigenvalues. 

With A. Melin [10]: we considered an /i-pseudodifferential operator with 
leading symbol p(x, ^) that is bounded and holomorphic in a tubular neigh- 
borhood of R^ in C* = C^ X C|. Assume that 

R^ n p~^(0) ^ is connected. (3.1) 

On R^ we have b(x,OI > 1/C, for |(a;,OI > C, (3.2) 

for some C > 0, 

(i3f?p(a:, ^), d'^p{x,^) are linearly independent for all (x, ^) G p^^(O) fl R^. 

(3.3) 
(Here the boundedness assumption near cxd and (13. 2p can be replaced by a 
suitable ellipticity assumption.) It follows that p~^(0) fl R'^ is a compact 
(2-dimensional) surface. 
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Also assume that 

\{^p, Qp}\ is sufficiently small on p-^(O) n R^ (3.4) 

Here "sufficiently small" refers to some positive bound that can be defined 
whenever the the other conditions are satisfied uniformly. 

When the Poisson bracket vanishes on p~^{0), this set becomes a La- 
grangian torus, and more generally it is a torus. The following is a com- 
plex version of the KAM theorem without small divisors (cf T.W. Cherry 
[3](1928), J. Moser ^(1958)), 

Theorem 3.1 (IJO^ ) There exists a smooth 2-dimensional torus T C p~^(0)n 
C^ , close to p~^(0) n R^ such that a\^ = and Ij(T) G R, j = 1,2. Here 
Jj(r) := J X'dx are the actions along the two fundamental cycles 71, 72 C T , 

and a = "^-^ dC,j A dxj is the complex symplectic (2,0)-form. 

Replacing p hj p — z for z in a neighborhood of G C, we get tori 
T{z) depending smoothly on z and a corresponding smooth action function 
I{z) = (/i(r(z)), /2(r(z))), which are important in the Bohr- Sommerf eld rule 
for the eigen-values near in the semi-classical limit h —>■ 0: 

Theorem 3.2 (IJm ) Under the above assumptions, there exists Oq G (|Z)^ 
and e{z; h) ^6^ + ei{z)h + e2{z)h^ + .. zn C°°(neigh (0, C)), such that for z 
in an h-independent neighborhood of and for h > sufficiently small, we 
have that z is an eigenvalue of P = p{x, hD^) iff 

li^ = k- e{z; h), for some k G Z^. (BS) 

Recently, a similar result was obtained by S. Graffi, C Villegas Bas [T3] . 

An application of this result is that we get all resonances (scattering 
poles) in a fixed neighborhood of G C for —h^A + V{x) if V is an analytic 
real potential on R^ with a nondegenerate saddle point at x = 0, satisfying 
V{0) = and having {{x,C,) = (0,0)} as its classically trapped set in the 
energy surface {p{x,^) = 0}. 



3.3 Diophantine case 

In this and the next subsection we describe a result from [SU] and the main 
result of [29] about individual eigenvalues for small perturbations of a self- 
adjoint operator with a completely integrable leading symbol. We start with 
the case when only Diophantine tori play a role. 
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Let P^{x,hD]h) on R^ have the leading symbol Pe(x,^) = p{x,^) + 
ieq{x, ^) where p, q are real and extend to bounded holomorphic functions on 
a tubular neighborhood of R^. Assume that p fulfills the ellipticity condition 
(13.21) near infinity and that 

P,=o = Pix,hD) (3.5) 

is self-adjoint. (The conditions near infinity can be modified and we can also 
replace R^ by a compact 2-dimensional analytic manifold.) 

Also, assume that Pe{x, ^; h) depends smoothly on < e < eo with values 
in the space of bounded holomorphic functions in a tubular neighborhood of 
R^, and P^ ~ Pe + ^Pi,e + h'^P2,e + ■■■, when /i — > 0. 

Assume 

p~^{0) is connected and dp ^ on that set. (3.6) 

Assume complete integrability for p: There exists an analytic real valued 
function / on T*R^ such that Hpf = 0, with the differentials df and dp being 
linearly independent almost everywhere on p~^(0). {Hp = p'^ ■ -^ — p'^ ■ -^ is 
the Hamilton field.) 

Then we have a disjoint union decomposition 

p-\0)nT*R^ =[JA, (3.7) 

AeJ 

where A are compact connected sets, invariant under the Hp flow. We as- 
sume (for simplicity) that J has a natural structure of a graph whose edges 
correspond to families of regular leaves; Lagrangian tori (by the Arnold- 
Mineur-Liouville theorem). The union of edges J\S possesses a natural real 
analytic structure. 

Each torus A G J \ S* carries real analytic coordinates Xi,X2 identifying 
A with T^ = RV27rZ2, so that along A, we have 

H. = a.^^+a,£-, (3.8) 

where ai,a2 G R. The rotation number is defined as the ratio uj{A) = [ai : 
02] e RP\ and it depends analytically on A G J \ 5*. We assume that ci;(A) 
is not identically constant on any open edge. 

We say that A G J \ S" is respectively rational, irrational, diophantine if 
ai/a2 has the corresponding property. Diophantine means that there exist 
a > 0, c? > such that 

\{a,,a2)-k\>--—^, 0^keZ\ (3.9) 
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We introduce 

{q)T = 7^ / go eMtHp)dt, T>0, (3.10) 

^ J-T/2 

and consider the compact intervals Qoo{A.) C R, A G J, defined by, 

<5oo(A) = [ lim inf(g)T, lim sup(g)r]. (3.11) 

T— »oo A T— >cx) yY 

A first localization of the spectrum a{Pe{x, hD^; h)) ([30]) is given by 

53((t(P,) n {z; \^z\ < 6}) C e[inf |J Q^{A) - o(l),sup |J Qoo(A) + o(l)], 

AeJ AeJ 

(3.12) 

when 6,e,h -^ 0. 

For each torus A E J\S, we let (g)(A) be the average of g|^ with respect 
to the natural smooth measure on A, and assume that the analytic function 
J \ S 3 A \-^ {q){A) is not identically constant on any open edge. 

By combining (13. 8p with the Fourier series representation of q, we see that 
when A is irrational then Qoo(A) = {(g) (A)}, while in the rational case, 

Qoo(A) C (g)(A) + 0(^^^^^)[-l, 1], (3.13) 

when ci;(A) = — and m G Z, n G N are relatively prime. 

Let Fo G UAejQoo(A) and assume that there exists a Diophantine torus 
A^ (or finitely many), such that 

(g)(Arf) = Fo, rfA(g)(A,)^0. (3.14) 

With M. Hitrik and S. Vd Ngoc we obtained: 

Theorem 3.3 (f3^) Assume also that Fq does not belong to Qoo{A) for 
any other A G J. Let < 6 < K < oo. Then 3C > such that for 
h> small enough, and k^ < e < h^ , the eigenvalues of P^ in the rectangle 
\'^z\ < h^/C, \^z - e?R:Fo\ < eh^/C are given by 

P^^\h{k -^)-f, e; h) + 0{h^), k G Z^ 
Here P^°°\^, e; h) is smooth, real-valued for e = and when h -^ we have 

oo 

P(°°)(e,6;/i) ~5^/i^p(°^)(e,e), p(-) =p^(0 + z6(g)(0 + O(e^), (3.15) 



corresponding to action angle coordinates 

19 



In [30] we also considered applications to small non-self-adjoint pertur- 
bations of the Laplacian on a surface of revolution. Thanks to (13.131) the 
total measure of the union of all Qoo(A) over the rational tori is finite and 
sometimes small, and we could then show that there are plenty of values Fq, 
fulfilling the assumptions in the theorem. 

3.4 The case with rational tori 

Let Fo be as in (13.141) but now also allow for the possibility that there is a 
rational torus (or finitely many) A^, such that 

Fo e Qoo(Ar), Fo ^ (g)(A,), (3.16) 

dA((g))(A.) ^ 0, dA{Lu)iAr) ^ 0. (3.17) 

Assume also that 

Fo ^ goo(A), for all A G J \ {A^, AJ. (3.18) 

With M. Hitrik we showed the following result: 

Theorem 3.4 (J2^) Let 6 > be small and assume that h <^ e < h^^^ , or 
that the subprincipal syvahol of P vanishes and that h^ <^ e < Ha^^ . Then 
the spectrum of P^ in the rectangle 

is the union of two sets: E^ U Ej., where the elements of E^ form a distorted 
lattice, given by the Bohr-Sommerfeld rule Ii3.15\) . with horizontal spacing 
X h and vertical spacing x eh. The number of elements i^{Er) of Er is 
0{e'l^/h^). 

NB that #(Ed) X ei+V/i2. 

This result can be applied to the damped wave equation on surfaces of 
revolution. 



3.5 OutUne of the proofs of Theorem 13.31 and 13.4 

The principal symbol of P^'m p^= p + ieq + (9(e^). Put 

rT/2 
-T/2 



1 f^/^ 

{q)T = 7f ? ° exp{tHp)dt. 

^ J-TI2 
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As in Section [2] we will use an averaging of the imaginary part of the symbol. 
Let J{t) be the piecewise affine function with support in [— |, |], solving 

J'(t) = 5(t)-l[_i,i](t), 

and introduce the weight 



t 

T' 



Grit) = I J{--)qoexp{tHp)dt. 



Then HpCx = <? — (q')T, implying 

p, o e^ieHc^) =p + te{q)T + Orie^). (3.19) 

The left hand side of f l3.19p is the principal symbol of the isospectral 
operator e~t'^T(a;,/iDa;)Qp^Qg|-GT(a;,/iDa;) g^j^^ under the assumptions of Theorem 

13.31 resp. 13.41 its imaginary part will not take the value ieFQ on p^^(O) away 
from Arf resp. A^UA^. This means that we have localized the spectral problem 
to a neighborhood of A^ resp. A^ U A^. 

Near A^ we choose action-angle coordinates so that A^ becomes the zero 
section in the cotangent space of the 2-torus, and 

Pe{x, = PiO + mix, + 0{e^). (3.20) 

We follow the quantized BirkhofF normal form procedure in the spirit of 
V.F. Lazutkin and Y. Colin de Verdiere [SSlll]: solve first 

i7pG = g(x,0-(g(-,0), (3.21) 

where the bracket indicates that we take the average over the torus with 
respect to x. Composing with the corresponding complex canonical trans- 
formation, we get the new conjugated symbol 

M0+^e(g(-,0) + O(e' + r). 

Here the Diophanticity condition is of course important. 
Iterating the procedure we get for every A^, 



p,oexp(ifc,^))=p(0+^e((g)(0 + (^(e,0)+(^((e,e 



,Af+l^ 



independent of x 



This procedure can be continued on the operator level, and up to a small error 
we see that P^ is microlocally equivalent to an operator P^{hD^, e; h). At least 
formally. Theorem 13.31 then follows by considering Fourier series expansions, 
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but in order to get a full proof we also have take into account that we have 
constructed complex canonical transformations that are quantized by Fourier 
integral operators with complex phase and study the action of these operators 
on suitable exponentially weighted spaces. 

Near A^ we can still use action-angle coordinates as in f l3.20p but the 
homological equation (13.211) is no longer solvable. Instead, we use secular 
perturbation theory (cf the book [2Z]), which amounts to making a partial 
Birkhoff reduction. 

After a linear change of x-variables, we may assume that p{^) = ^2 + 0{C,'^) 
and in order to fix the ideas = 1^2 + ^i- Then we can make the averaging 
procedure only in the a;2-direction and reduce p^ in (13.201) to 

Mx,0 = ^2 + ei+O{e)+O{{e,0n, 

'• V ' 

independent of X2 , 

»e2+«i+"(9>2(^l>€) 

where ((?)2(a^i,0 denotes the average with respect to X2- 

Carrying out the reduction on the operator level, we obtain up to small 
errors an operator Pe{xi, hD^^-, hD^^; h) and after passing to Fourier series in 
X2, a family of non- self-adjoint operators on S'^^: Pe{xi, hD^^, hk; h), A; G Z. 

The non-self-adjointness and the corresponding possible wild growth of 
the resolvent makes it hard to go all the way to study individual eigenval- 
ues. However, it can be shown that in the region |^i| 3> e^/^ (inside the 
energy surface p = 0) we can go further and (as near A^) get a sufficiently 
good elimination of the x-dependence. This leads to the conclusion that the 
contributions from a vicinity of A^ to the spectrum of P^ in the rectangle 

e e^+^ 

\^z\ < -, l^^-eFol < 



come from a neighborhood of A^ of phase space volume (9(e^/^). 

This explains heuristically why the rational torus will contribute with 
O^e^/"^ /h?) eigenvalues in the rectangle. 

The actual proof is more complicated. We use a Grushin problem reduc- 
tion in order to reduce the study near A^. to that of a square matrix of size 
0[e'l'^ lh?\ However, even if we avoid the eigenvalues of such a matrix, the 
inverse can only be bounded by 

expO(e3/V/i'). (3.22) 

What saves us is that away from A.,. U A^, we can conjugate the operator with 
exponential weights and show that the resolvent has an "off-diagonal decay" 
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like exp(— l/(C/i)). This implies that we can confine the growth in (13.221) to 
a small neighborhood of A^, if 

leading to the assumption e <C h'^^'^ in Theorem I3.4[ 



4 Weyl asymptotics for non-self-adjoint op- 
erators 

4.1 Introduction 

For self-adjoint differential (pseudo) differential operators we have (under 
suitable assumptions) the Weyl law for the asymptotic distribution of eigen- 
values, established in higher dimensions by H. Weyl [56j| in 1912 in the case 
of second order elliptic boundary value problems. 

In the semiclassical setting such results were obtained by J. Chazarain, 
B.Helffer-D. Robert, V.Ivrii and many others (see [TT] and further refer- 
ences there). Under suitable additional assumptions it states that if P = 
P^{x, hDx] h) is a self-adjoint /i-pseudodifferential operator with leading (real) 
symbol p(x, ^), then if fi C C is a domain intersecting R along a bounded 
interval, the number of eigenvalues of P in fi (intersected with R) satisfies 

#(a(P) n fi) = ^^(voi {p-'m) + 0(1)), h -. 0. 

A simple example is provided by the semiclassical harmonic oscillator 
^{{hD^Y+x^) on the real line which has the eigenvalues (A; + |)/i, A; = 0, 1, .... 

In the non-self-adjoint case Weyl-asymptotics is known to hold in some 
cases close to the self-adjoint case or for normal operators. 

We do not always have Weyl-asymptotics: Following Davies and Boul- 
ton (see [7]), we can consider the non-self-adjoint harmonic operator: P = 
^{{hD^y + ix^) whose eigenvalues are given by e*'^/'^(/c + \)h, k eN ([46j): 

The set of values of p = |(^^ + ix'^) is the closed first quadrant and if we 
choose the open bounded set Q to intersect the 1st quadrant but not the line 
'isz = ^z, we get vol {p~^{Q)) > 0, while there are no eigenvalues in Q. 

More generally, /^-differential operators with analytic coefficients often 
have their spectrum determined by complex-geometric quantities, and are 
likely not to obey the Weyl law. In particular, in the one dimensional case it 
often happens that the eigenvalues are concentrated to certain curves with 
branch points. 
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As we have seen in Theorem 11.31 we are often confronted with the pseu- 
dospectral phenomenon: On the image of p the resolvent may be very large 
even far from the spectrum. This causes the eigenvalues to be very sensitive 
to small perturbations of the operator (by Theorem 11.21) . 

In her thesis M. Hager (see [E]) considered a class of perturbed h- 
pseudodifferential operators on the real line of the form Ps = P{x, hD; h) + 
5q^{x), where P is analytic and q^^ is a random linear combination of the 
C /h first eigen-functions of an auxiliary operator. She showed that with 
probability very close to 1 when h ^ 0, Ps obeys Weyl asymptotics. 

Here, we shall discuss a generalization to the multidimensional case ob- 
tained with Hager [21]. The results will be much more general in many ways, 
but the class of perturbations will be slightly different. 

4.2 The result 

a) The unperturbed operator. Let m{p) > I, p = {x,C,) be an order function 
on R^", so that 

0<m(p)<Co(p-^)^M/"), 

we may assume that m G S{m) = {u e C°°(R2"); d'^u = Oaim), \/a G N^"}. 
Assume m>l and let 

P{p] h) ~ p{p) + hpi{p) + ... in S{m). 

Assume 3^0 G C, Co > such that 

\p{p) - ^ol > m{p)/Co (ellipticity). 

Let 



S=p(R2n)=p(R2«)uSoo, 

Soo = { hm p{pj); R^" 3 pj ^ oo} 



We write: P = P'"(x, hD;h). 

Let i7 CC C \ Soo be open and simply connected containing zq. Then 
using the pseudodifferential calculus, it is easy to show: 

1) cr(P) n i7 is discrete when h > is small enough. 

2) Ve > 0, 3h{e) > 0, such that ct(P) n fi C S + D{0, e),0<h< h{e). 

b) The random pertubation. Let < m, m < 1 be square integrable 
order functions, one of which is integrable. Let 5* G S{m), S G S{m) 
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be elliptic symbols. The corresponding operators are Hilbert-Schmidt with 
||^||hs,||^||hs = C?(/^-"/'). Let 

Qoj = S o^ aj,fc(c^)eje^ o S, 

where dj^k, j, A; G N, are independent complex Gaussian random variables 
with expectation value and variance 1, and (e,)f^ and (e,)f^ are orthonormal 
bases in L^(R"), CjClu = {u\ek)ej. 

Let M = Ci/i-" with C^ > 1. Then with probability > 1 - Ce''"'^"/^ 
the Hilbert-Schmidt and trace class norms of Q fulfil: 

||Q||hs<M, ||g||tr<Mi (4.1) 

Let r CC f2 be open with smooth boundary. 

Theorem 4.1 (fST]). Assume 

p{p) E dr ^ dp{p), dp{p) are linearly independent. (4-2) 

Let e, 6 depend on h with < e ^ 1, 

(implying that e > Const. /i In j). Then with probability > 1 — -^e W^^W^ ^ we 
have 

m-(Ps) n r) - ^voi {p-\m < c^.. (4.3) 

There is a similar result giving (14.31) simultaneously for all F in a suitable 
family. 

The assumption (14. 2 p implies that 9rn(9E = 0, so if we want to count the 
eigenvalues near 9S, we need to weaken that assumption. For z G neigh (dT), 
we put 

Kit) = vol {p G R2"; |p(p) - z|2 < t}. (4.4) 

Introduce the assumption 

3k g]0, 1], such that V,{t) = C(t'"), uniformly for z G neigh (dV), < t < 1. 

(4.5) 
Example. (S2D => (USD with k = 1. 
Example. The best that can happen when dT fl 9S 7^ is that 

p{p) edV^ {p,p}ip) ^ or {p, {p,p}}ip) ^ 0. (4.6) 

It is easy to see that (14.61) implies (14.51) with k = 3/4. (14.61) holds for the 
non-self-adjoint harmonic oscillator when ^ dT. 
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Theorem 4.2 (l21^) We assume U-^ - Let e, 5 depend on h with < e ^ 1, 
implying that e > Const. /i*^ In ^. Then for < r <^ 1 we have with probability 

— r ' 



|#(ffW)nr)----;voi(p-i(r))|< (4,7) 






+ Cn [r^ + ln( J)vol {p"\dV + D(0, r)))\\ 



for every fixed A^ e N . 

If K > |, we have Mo\{p-^{dV + £'(0,r))) = C(r2«-i) with 2^; - 1 > 
and in all cases we may assume that ln(i)vol (p~^((9r + D(0,r))) = (9(r"''), 
where ao > 0. Then we choose iV ^ 1, r = ei/(^+"o) and the right hand side 
of dHH) becomes C»(l)e"o/(i+"o)/i-'^. 

Again we have a similar theorem where the conclusion (14 .Zp is valid si- 
multaneously for all r in a suitable family. 

Recently, the author obtained similar results when Q^ is an operator of 
multiplication, see 



4.3 Outline of the proofs 

We can construct P G S{m) such that P(p; h) = P{p] h) for \p\ ^ 1 and 
\P{p;h) — z\ > m{p)/C for p G R^", z ^ Q. The eigenvalues of P in f2 
coincide with the zeros of the holomorphic function 

F{z;h) = detP,, (4.8) 

P, = {P{x,hD;h)-zy\P{x,hD;h)- z) 

The same remark holds for Ps and Fs defined as in (14.81) with P replaced by 
Ps, provided that (gj]) holds. 

For z G neigh ((9r), put Q = P*Pz- Let la{E) = max{E,a), where 
a = Ch, C ^ 1. Using semiclassical analysis, we can show that under the 
assumption igB) (cf p]) 

lndetQ<lndetl„(Q) = , \, { ( ( \nqdxdi + 0{l)h''\n\), (4.9) 

(27rft,)" J J h 

where q = \pz\'^ is the leading symbol of Q. 
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Since 

In det g = In det P*P, = In \F{z; h)\'^ 

we have 

\n\F{z;h)\<j^{JJ\n\p,\dxd^ + 0{l)h^\n^), (4.10) 

For S > small enough, we get the same upper bound for In \Fs{z; h)\ (pro- 
vided that (O) holds). 

The main step in the proof is to get a corresponding lower bound for 
each fixed z with a probability close to 1. In the multidimensional case this 
boils down to a question about random determinants. Let z G neigh (SF). 
Let ei, 62, ... be the first eigen-functions oi Q = P*Pz and let /i, /2, ... be the 
first eigen-functions for PzPz- The two operators have the same eigenvalues 
< Ai < A2 < ... 

We can arrange so that 



Let A^ = N{a) = 4^{j',\j < a} {a = Ch, C ^ 1). Semiclassical analysis 
gives that 

N = C(/i'^""). 

Consider 



r^= { ^^ -^ I -.L'xa' -^L'xa\ 



N 



1 



V^ has an inverse 



with El,E^_,E'^^ "explicit", and 

In I detPY = A^lna + det laiP^Pz)- 
For Ps = P + SQ^, we form 

Pf = {P- z)-\P ^z + 6Q^) = Pz + 6Q, 
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and since 5||(5t^|| < SCoh " <^ 1, 

-p<5 .. 

is invertible with inverse 

£'=("::. r/ 1 ^ s'>. 

Here 

^i+ = ^°+ + 5E^_Q^El + "small", (4.11) 

and we can show by perturbative arguments that 

In det V^ = In det P° + 0(^M^^/^), 

Ja 



/ 


Pt 


i?_ 


v 


R+ 





E' 


E'. 





leading to 

ln|detP^| = .^ \, ( / / ln|p,|rfxrfe + C(/i''ln^)). (4.12 



On the other hand, computations in fSOJ can be used to get 

In I det Pf I = In I det P^ I + In I det Ei^ \ (4. 13) 

Using (14.111) . we can view E^_^_ as a random matrix of size 0{h'^~'"'), close in 
a suitable sense to one with independent Gaussian random variables as its 
entries. This can be used to show: 

For every z G neigh (SF), we have a nice lower bound on In | det i?i_,_| 
with probability close to 1. (14.131) then gives a corresponding lower bound 
on In I det P^ I- 

To complete the proof of Theorem 14.11 we can apply the following result 
of M. Hager [HI [18] with h = h"^: 

Proposition 4.3 ([Tl[ [7g|/). Let F and fi be as above. Let (f) G C{^; R) be 
smooth near dT. Let f = f{z; h) be holomorphic in Q with 

\f{z;h)\ < e'^(")/'^, z G neigh ((9F), < /J < 1. 

Assume there exist e = e{h) <^ 1, Zk = Zk{h) & Q, k E J = J{h), such that 

dTG[JD{z,,^e), #J = 0(±), 
ke.J ^^ 

|/(^;t;/i)|>e(*(^'=)-^)/^ A;G J. 
28 



Then, 

#(/-i(o) n r) = i [ [ iA<i^)d{^z)d{Qz) + o{^). 

2nh J Jr h 

For the proof of Theorem 1.2 we use an improved version of this result, 
see [TIJ. 

4.4 Comparison with Theorem 13.21 

From the example with the non- self-adjoint harmonic oscillator in dimension 
1, we have seen that Weyl asymptotics does not always hold for differential 
operators in one dimension, when the coefficients are analytic. If we add a 
small random perturbation to the non-self-adjoint harmonic oscillator, the 
theorems above and the main result in [IS] show that with probability close 
to 1 the eigenvalues will no longer be confined to a half-line but will tend to 
fill up the range of the principal symbol p with a density that is given by 
{27r h)~"'p^{dv), where dv denotes the symplectic volume element on R^" and 
p*{dv) is the direct image under p. 

From this simple one-dimensional example it is easy to build examples in 
higher dimension when Weyl asymptotics does not hold. In the 2-dimensional 
case, we can also consider the situation when the unperturbed operator P 
satisfies the assumptions of Theorem 13. 2[ It is then natural to compare the 
distribution law given by Theorem 13.21 for P and the one given by the the- 
orems 14.11 14.21 for the random perturbations. To leading order in h, we get 
Weyl asymptotics already for P in the (close to normal) case when {p,p} 
vanishes identically. In general however, we get different asymptotic distri- 
butions already to leading order ([^8]). 
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